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ABSTRACT 

Based on the detected variability time scales of X-ray and TeV gamma-ray emis- 
sion, and the observed multiwavelength photon spectrum, of Mrk 421 we place con- 
straints on the allowed parameter space (magnetic field and Doppler factor of the 
emission region) for the homogeneous synchrotron self-Compton model. The spectra 
calculated for the allowed parameters are marginally consistent with the available spec- 
tral information above ~ 1 TeV reported by the Whipple Observatory in the case of 
a 1 day flare time scale. However, for the recently reported very short duration flares 
varying on a time scale of 15 min, the calculated spectra are significantly steeper, 
suggesting that the homogeneous synchrotron self Compton model has problems in 
describing the relatively flat observed spectra extending above a few TeV. We deter- 
mine the maximum ratio of TeV gamma-ray luminosity to X-ray luminosity during 
flaring which is allowed by the homogeneous synchrotron self-Compton model for the 
case of no significant photon-photon absorption in the source. 
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1 INTRODUCTION 

Very high energy (VHE) 7-ray emission has been detected 
in recent years from two BL Lac objects Mrk 421 (Punch 
et al. 1992, Petry et al. 1996) and Mrk 501 (Quinn et 
al. 1996). The emission varies significantly on different time 
scales, from weeks and days (Kerrick et al. 1995, Schubnell 
et al. 1996, Buckley et al. 1996), to fractions of an hour 
(Gaidos et al. 1996). Multiwavelength observations of Mrk 
421 show that the TeV 7-ray flares are simultaneous with 
the X-ray flares observed by the ASCA satellite, and that 
the power emitted in these two energy ranges is compara- 
ble (Takahashi et al. 1996, Macomb et al. 1995, Buckley et 
al. 1996). However during this same period, the amplitude 
of variations in the lower energy emission (UV-optical) and 
the 7-ray emission in the EGRET energy range was much 
less than that of the very high energy 7-ray emission (Ma- 
comb et al. 1995, Buckley et al. 1996, Lin et al. 1994). VHE 
emission from Mrk 421 has been observed up to ~ 4 TeV in 
the quiescent state (Mohanty et al. 1993) , and up to at least 
~ 8 TeV in the high state (Krennrich et al. 1997). Emis- 
sion extending to similar energies has also been observed 
recently from Mrk 501 (Breslin et al. 1997). Summing the 
HEGRA data from several nearby blazars (including Mrk 
421) there is some evidence (6.5<r) of 7-ray emission above 



50 TeV (Meyer & Westerhoff 1996, Meyer 1997), and for 
Mrk 421 alone the significance is 3.8<r. 

Gamma ray emission from active galactic nuclei (AGN) 
is often interpreted in terms of the homogeneous "syn- 
chrotron self-Compton model" (SSC) in which the low en- 
ergy emission (from radio to X-rays) is synchrotron radiation 
produced by electrons which also up-scatter these low energy 
photons into high energy 7-rays by inverse Compton scat- 
tering (ICS) (Macomb et al. 1995, Inoue & Takahara 1996, 
Bloom & Marscher 1996, Mastichiadis & Kirk 1997). In this 
model all the radiation comes from this same region in the 
jet. Such a picture can naturally explain synchronized vari- 
ability at different photon energies. More complicated (inho- 
mogeneous) SSC models are also proposed which postulate 
that the radiation at different energies is produced in differ- 
ent regions of the jet (e.g. Ghisellini et al. 1985, Maraschi 
et al. 1992). It has also been argued that the 7-ray emission 
from Mrk 421 can be explained by electrons scattering in 
the Klein-Nishina regime (Zdziarski & Krolik 1993). 

Photon-photon pair production on the infrared back- 
ground radiation was expected to prevent observation above 
~ 1 TeV. The observation of 7-rays from Mrk 421 up to 
8 TeV without a cut-off in the spectrum, and certainly up to 
50 TeV, was unexpected based on calculations for reasonable 
models of the infrared radiation field (see e.g. Stecker and 
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De Jager 1997 and references therein). Even if cascading is 
included in the infrared background (Protheroe and Stanev 
1993, Entel and Protheroe 1995) observation of 7-rays of 
these energies seemed unlikely. However, the infrared back- 
ground is not well known and it may be possible to observe 
the nearest blazars up to energies somewhat below ~ 100 
TeV where absorption on the cosmic microwave background 
will give a sharp cut-off. In our calculations below, we shall 
therefore neglect photon-photon pair production on the in- 
frared background. 

The purpose of this paper is to confront the homoge- 
neous SSC model with the results of recent observations and, 
if possible, to derive the parameters of the emission region in 
the jet from which this radiation originates, i.e. its Doppler 
factor and magnetic field strength. 



2 CONSTRAINTS ON A HOMOGENEOUS SSC 
MODEL 

Let us consider relativistic electrons confined in a "blob" 
which moves along the jet with the Doppler factor D and 
has magnetic field B. In the homogeneous SSC model the 
radii of the emission regions of low energy photons (r;), X- 
ray photons (rj), and TeV 7-rays (r 7 ) are the same. This 
region is constrained by the variability time scale observed, 
e.g. in TeV 7-rays, t var (s), 



n = r 1 = r x 



0.5cDt v 



(1) 



The differential photon density in the blob frame of low en- 
ergy synchrotron photons (phot. MeV -1 cm -3 ) is then given 
by 



i(e') 



4d 2 F{e) 

r 3+2 734 ■ 



(2) 



where d w 187 Mpc is the distance to Mrk 421 (for H = 50 
km s _1 Mpc -1 , and z — 0.031), e = De' and e are the 
photon energies in the observer's and the blob rest frames, 
and c is the velocity of light. 

The differential photon flux in the optical to X-ray re- 
gion observed from Mrk 421 during the 16 May 1994 flare 
(photons cm -2 s _1 MeV -1 ) can be approximated by a bro- 
ken power-law, 



F(e) 



he 
hie 
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10~ 3 e 6 <e<e b , 



(3) 



where /3i = 1.8, #2 = 2.3, 61 = 6.9 x 10~ 4 and b 2 = 2.8 x 
10~ 5 are obtained from ASCA observations, and e b = 1.65 x 
10" 3 MeV is the energy at which a break in the synchrotron 
spectrum is observed. Eq. ^ is based on the peak 2-10 keV 
luminosity and spectral index given in Fig. 3 of Takahashi 
et al. (1996) (for e > e b ), and from Fig. 2 of Macomb et al. 
(1995) we estimate the spectrum for e < e b . Note that for at 
least 2 decades below 10 _3 ej, the spectrum is uncertain. 

The shape of synchrotron spectrum defines the shape 
of the electron spectrum in the blob rest frame which can 
be approximated by 



dN 

drf 



a l7 '- a i 0.032 7 £ < 7' < 7 i, 



where 7' is the Lorentz factor in the blob frame, 
75 = (2e b /De B ) 1/2 , 



(4) 



(5) 



qi = 2.6, a 2 = 3.6, e B = m e c 2 B/B cr , B cr = 4.414 x 10 13 
G, m e is the electron rest mass, ai = 02/7^, and 02 can 
be obtained from fitting the observations. Note that below 
10~ 3 ^ 2 -f' b = 0.0327^ the spectrum is uncertain. 

The spectrum of Mrk 421 shows two clear bumps which, 
during the outburst stage, extend up to at least ~ 10 keV 
(Takahashi et al. 1996), and ~ 8 TeV (Krennrich et al. 1997). 
These multiwavelength observations of Mrk 421 allow us to 
define the ratio 77 of the power emitted at a 7-ray energy, 
E-f, at which the emission is due to Compton scattering, to 
the power emitted at an energy, e, at which the emission is 
due to X-ray synchrotron radiation, 



dedt 



dN ,2 dN 
'E'dt' ~ l>/( de'dt' 



(6) 



where the primed quantities are measured in the blob frame. 
For the power at 7-ray energies we adopt the value reported 



for the threshold of the Whipple telescope at E y 



0.3 



TeV (Macomb et al. 1996), and for the power at X-ray syn- 
chrotron energies we take the value corresponding to the 
peak emission at e = e b (Takahashi et al. 1996). For these 
two energies 77 ~ 1.2. 

The synchrotron spectrum at e' in the above formula 
(Eq. ^) can be obtained approximately analytically from the 
relation 



, dN 
de'dt 



de' 



dN 



jdy'bsyni-y'), 
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where dN/d'y' is the electron spectrum (Eq. j4J). The charac- 
teristic energy of synchrotron photons is given by 



e 0.5es7 ' 



(8) 



the energy loss rate of electrons is 6 syn (7 / ) = kl/B'f' 2 , where 
k = 4cctt/3, ot is the Thomson cross section, and Ub ~ 
2.5 x 10 4 _B 2 (MeV cm -3 ) is the magnetic field energy density. 
The synchrotron spectrum emitted by electrons with power- 
law spectral index a, multiplied by the square of the photon 
energy, is given by 



dN 
de'dt' 
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The ICS part of the Eq. ^ cannot be obtained analytically 
in the general case because of the complicated form of the 
Klein-Nishina cross section, and so we compute this numer- 
ically using 



F /2 _ ,2 
dE'dt' ' ' 



dN 



dN(j',E' 



dt'de'dE: 



^de'di, (10) 



where 7 ^ in « E'Jm e c 2 , e' min = fy/Wtf - E'^/m.c 2 )], 
El, = Ey/D, and dN (7' , E^) / dt' de' dE\ is the ICS spectrum 
(see Eq. 2.48, in Blumenthal & Gould 1970) produced by 
electrons with Lorentz factor 7' which scatter synchrotron 
photons in the blob having the spectrum given by Eqs. ^ 
and|[ 

Having determined the spectra in Eq. ^, we can now 
investigate the parameter space (magnetic field strength in 
the blob, B, and Doppler factor, D) for the homogeneous 
SSC model which is consistent with the value of r\ = 1.2. 
In Figs. |l|(a) and [j](b) we show the allowed value of B as a 
function of D (thick full curves) for the case of outbursts as 
reported by the Whipple Observatory which varied on (a) a 
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- 1 day (Buckley et al. 1996, Schubnell et al. 1996) and (b) 
a ~ 15 min (Gaidos et al. 1996) time scale. 

2.1 Limits from variability time scales 

The variability time scales observed in TeV 7-rays by the 
Whipple Observatory, and the reports of simultaneous flares 
observed in X-rays by ASCA allow us to place a further 
constraint on the homogeneous SSC model. A significant 
decrease in the observed TeV 7-ray and X-ray fluxes may 
only occur if the electrons have sufficient time to cool during 
the flare, 

t cool < t vaI D. (11) 

The cooling time scale for synchrotron losses of electrons 
with Lorentz factor 7^, which contribute mainly to syn- 
chrotron photons at the peak of the spectrum, is given by 



,/syn 
^ cool 



kUsi h ' 



(12) 



Eqs. (^), (|ll[), and (Q) allow us to place a lower limit on 
the magnetic field in the blob 



B > 15.lt 



-2/3-1/3 



D 



-1/3 



(13) 



We next estimate the ICS cooling time and require this 
to be less than the variability time scale. For the soft photon 
spectrum we adopt, some interactions with energetic elec- 
trons will be in the Thomson regime, and others will be in 
the Klein-Nishina regime with relatively small energy loss. 
To calculate the cooling time, we neglect interactions in the 
Klein-Nishina regime, i.e. with photons above e' T ~ m e c 2 /j' , 
and use the simple Thomson energy loss formula 

(14) 



/ICS 
v cool 



kU rad {< 4)7" 



where 



Urad{< e T ) 



n(e')e'de ~ 5Ce^!' 



(15) 



C = 4d 2 /(c 3 ^ ar D 4+ ' 31 ), and we have used the fact that e' T < 
e' b for electrons emitting at TeV energies. 

Since the recent report on the Mrk 421 flare observa- 
tions by the Whipple Observatory gives no evidence of a 
spectral break at high energies, we assume that the break in 
the 7-ray spectrum is below but close to the threshold of the 
Whipple observations, i.e. at ~ 0.3 TeV. This in turn means 
that electrons with Lorentz factor 7^ must cool during the 
7-ray flare, and so we shall use 7' = 7^ in Eq. ([[4]) to esti- 
mate the ICS cooling time scale. Then, from Eqs. (|ll]) and 
( p^ ) we can obtain an upper limit on the magnetic field in 
the blob as a function of Doppler factor, as required by the 
ICS cooling argument 

(16) 
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Figure 1. 

The parameter space (B, D) allowed by the homogeneous SSC 
model for variability in Mrk 421 on a time scale t V ar of (a) 1 day 
and (b) 15 min. The thick full curves show the condition for r\ = 
1.2, and the thin full curves (labelled by value of rf) show the 
condition for other values of r\. The other curves give allowed 
ranges for: efficient cooling of electrons during the flare time scale 
by synchrotron radiation (Eq. |li| ) - area above dot-dash curve 
labelled 'Synch'; efficient cooling by ICS (Eq. |li| ) — area below 
dotted curve labelled 'ICS'; escape of 8 TeV 7-rays — area to right 
of long-dashed curve; escape of 50 TeV 7-rays - area to right of 
short-dashed curve. The shaded area is the allowed region for 
the parameters for a spectrum extending to 8 TeV. The marginal 
values of (B, D) which just fulfil the condition r) = 1.2 are marked 
by squares and labelled (i) to (iv). 



2.2 Absorption of gamma-rays in the blob 
radiation 



The constraints on B and D of the blob, derived above 
(Eqs. and ^|), are shown in Figs. |l](a) and [j](b) for the 
two variability time scales. In the next section, we derive an 
additional limit on the Doppler factor of the blob which can 
be obtained based on the non-observation of absorption of 
7-rays by photon-photon pair production with soft photons 
in the blob radiation. 



The observation of 7-ray flares with a spectrum extending up 
to ~ 8 TeV, or even 50 TeV, allows us to place a lower limit 
on the Doppler factor of the blob under the assumptions 
of the homogeneous SSC model. Using the observed soft 
photon spectrum of Mrk 421 (Eq. ^|) we can compute the 
optical depth t(E'j,D) for 7-ray photons with energy 
for e pair production inside the blob, 
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Figure 2. 

The adopted soft photon spectrum (solid curve at 10 -9 to 
10~ 4 GeV), and the ICS spectra calculated in terms of the ho- 
mogeneous SSC model for the four marginal values of (B, D) 
of the allowed regions marked by the numbered open squares in 
Figs. p|(a) and pj(b): (i) solid curve; (ii) dotted curve; (iii) short- 
dashed curve; (iv) long-dashed curve. Calculated spectra are com- 
pared with the observations of Mrk 421 (Macomb et al. 1995) 
during flaring (asterisk) and in a quiscent state (open squares). 



r(E'D) 



de"- 



dssa(s), (17) 



where n(e') is the differential photon number density and 
cr(s) is the total cross section for photon-photon pair produc- 
tion (Jauch and Rohrlich 1955) for a centre of momentum 
frame energy squared given by 



s = 2e'E!y{l - cos 6») 



(18) 



where 9 is the angle between the directions of the energetic 
photon and soft photon, and 



Smin — (2/n e c ) , 
(2m e c 2 ) 2 
£min - 4EI, : 

The condition 



(19) 
(20) 
(21) 

(22) 



t(b;,d) < i., 

gives us the limit on D. These lower limits are shown by the 
dashed lines in Figs. |l|(a) and [j](b) for variability time scales 
ivar = 1 day and 15 min, and for _E 7 = 8 TeV (long-dashed 
lines, Krennrich et al. 1997) and 50 TeV (short-dashed lines, 
Meyer & Westerhoff 1996). The allowed regions of the pa- 
rameter space {B, D) determined by Eqs. [H| and ^ 
are shown in Figs. [j](a) and [j](b) by the shaded areas. 



principle produce flares with r\ = 1.2 as required. We note 
that the values of (B, D) used in earlier modeling of the 
Mrk 421 spectrum (Inoue & Takahara 1996, Mastichiadis & 
Kirk 1997, Stecker, De Jager & Salamon 1996) are generally 
consistent with the parameter space derived by us. In order 
to determine if the broad band spectrum expected in the 
homogeneous SSC model is consistent with the 7-ray ob- 
servations during flaring, we compute the synchrotron and 
ICS spectra for four example parameters (B, D) from the 
allowed region indicated by points (i) to (iv) in Figs. |l](a) 
and 0(b). The calculated spectra are shown in Fig g. Note 
that in each case, the lowest energy we predict corresponds 
to 7' = O.O3275 (see Eq. ^| and comments below) which 
depends on the magnetic field. For the 4 cases this gives 
E' y « 36 GeV (i), 120 GeV (ii), 19 GeV (iii), and 47 GeV 
(iv), for the minimum energies for which we can predict the 
the 7-ray spectrum with any confidence in the homogeneous 
SSC model. 

For energies between 0.8 TeV and 8 TeV corresponding 
to observations made by the Whipple observatory during re- 
cent flaring of Mrk 421 (Krennrich et al. 1997), our predic- 
tions of the spectral index in the homogeneous SSC model 
range from 2.65 for 1 day variability and case (i), to 2.85 
for 15 minute variability and case (iv). The results obtained 
during flaring by Krennrich et al. (1997) up to ~ 8 TeV 
are consistent with the spectrum of Mohanty et al. (1993) 
taken during a quiescent state where the spectral index was 
~ 2.25 ± 0.19 ± 0.3 between 0.4 - 4 TeV. Given the error 
bars, this is just consistent with the spectral index of 2.65 
predicted for 1 day variability and case (i). However, we note 
that the calculated spectrum shows a break close to ~ 1 TeV 
which should be seen in the Whipple observations. 

In the case of a flare varying on a 15 min time scale, 
it seems that the spectra obtained in terms of the homoge- 
neous SSC model are not consistent with the relatively flat 
spectrum of the Whipple observations. The lower sensitivity 
HEGRA Cherenkov observations report a very steep spec- 
trum above ~ 1 TeV (spectral index 3.6 ± 1.) during the 
Dec. 94 - May 95 monitoring (Petry et al. 1996). However 
these observations refer not to outburst emission, but rather 
to quiescent emission since the spectrum is integrated over 
a long period. 

In conclusion, detailed spectral measurements in the en- 
ergy range above 0.3 TeV combined with the observations 
in the optical-X-ray range should allow one to determine 
precisely the parameters of the emission region (relativistic 
blob) and in general answer the question of the applicabil- 
ity of the homogeneous SSC model for 7-ray production in 
blazars. We note also that the absorption and synchrotron 
cooling conditions do not allow flares with 1 day variability 
having 77 > 200 (8 TeV) or 77 > 60 (50 TeV) - see the thin 
solid curves in Fig. |l](a). Similarly, for 15 minute variability 
T) > 70 (8 TeV) or 77 > 20 (50 TeV) are not allowed (see 
Fig. [lja). Observation of such huge 7-ray outbursts without 
accompanying X-ray outbursts would be inconsistent with 
the homogeneous SSC model. 



3 DISCUSSION AND CONCLUSION 

Inspection of the Figs. [j](a) and [l](b) shows that for some 
values of (B, D), i.e. the region of the thick full line in- 
side the shaded area, the homogeneous SSC model can in 
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